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ABSTRACT 

A  near-wall  two-equation  turbulence  model  of  the  K  —  e  type  is  developed  for  the  descrip¬ 
tion  of  high-speed  compressible  flows.  The  Favre- averaged  equations  of  motion  are  solved  in 
conjunction  with  modeled  transport  equations  for  the  turbulent  kinetic  energy  and  solenoidal 
dissipation  wherein  a  variable  density  extension  of  the  asymptotically  consistent  near-wall 
model  of  So  and  co-workers  is  supplemented  with  new  dilatational  models.  The  resulting 
compressible  two-equation  model  is  tested  in  the  supersonic  flat  plate  boundary  layer  -  with 
an  adiabatic  wall  and  with  wall  cooling  -  for  Mach  numbers  as  large  as  10.  Direct  compar¬ 
isons  of  the  predictions  of  the  new  model  with  raw  experimental  data  and  with  results  from 
the  K  —  uj  model  indicate  that  it  performs  well  for  a  wide  range  of  Mach  numbers.  The 
surprising  finding  is  that  the  Morkovin  hypothesis,  where  turbulent  dilatational  terms  are 
neglected,  works  well  at  high  Mach  numbers  provided  that  the  near  wall  model  is  asymptoti¬ 
cally  consistent.  Instances  where  the  model  predictions  deviate  from  the  experiments  appear 
to  be  attributable  to  the  assumption  of  constant  turbulent  Prandtl  number  -  a  deficiency 
that  will  be  addressed  in  a  future  paper. 
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tract  No.  NASI- 18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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1.  INTRODUCTION 


The  direct  numerical  simulation  of  compressible  turbulent  flows  -  at  high  Reynolds  num¬ 
bers  with  all  scales  resolved  -  will  not  be  possible  for  the  foreseeable  future,  if  ever  at  all. 
Turbulence  modeling  will  continue  to  play  a  crucial  role  in  the  computation  of  high-speed 
aerodynamic  flows  associated  with  the  design  of  advanced  aircraft.  In  these  applications, 
near-wall  turbulence  modeling  is  extremely  important  for  the  accurate  prediction  of  wall 
transport  properties,  such  as  the  skin  friction  and  heat  transfer  coefficients,  which  are  pivotal 
for  design.  Despite  its  technological  importance,  progress  in  near-wall  turbulence  modeling 
has  been  slow.  Many  of  the  commonly  used  near-wall  models  -  which  typically  contain  a 
variety  of  ad  hoc  wall  damping  functions  -  are  not  asymptotically  consistent  and  yield  poor 
predictions  even  in  simple  incompressible  boundary  layers  (see  Patel,  Rodi  and  Scheuerer^, 
Myong  and  Kasagi^  and  Speziale,  .4bid  and  Anderson^).  These  deficiencies  can  be  fatal 
when  turbulence  models  are  applied  to  separated  flows,  or  other  complex  boundary  layers, 
that  require  the  governing  equations  to  be  integrated  directly  to  the  wall. 

High-speed  compressible  flows  present  a  whole  range  of  new  problems  to  near-wall  turbu¬ 
lence  modeling.  Shock/boundary  layer  interactions  with  turbulence  amplification  and  flow 
separation  represent  but  two  examples.  Two-dimensional  equilibrium  turbulent  boundary 
layers  for  supersonic  flows  are  less  of  a  problem  provided  that  the  external  Mach  number  Moq 
is  not  too  large.  For  these  flows  it  is  generally  believed  that  Morkovin’s  hypothesis'*  -  com¬ 
monly  interpreted  to  mean  that  the  turbulence  statistics  are  only  altered  by  compressibility 
effects  through  changes  in  the  mean  density  -  is  valid  for  at  least  the  range  0  <  A/oo  <  5. 
This  hypothesis  allows  for  the  use  of  variable  density  extensions  of  existing  incompressible 
turbulence  models  for  which  dilatational  effects  are  neglected.  However,  even  the  ability  of 
these  models  to  reliably  predict  mean  velocity  profiles  in  two-dimensional  equilibrium  bound¬ 
ary  layers  for  .Mach  numbers  M^c.  >  5  has  been  recently  called  into  question  (see  Bradshaw, 
Launder  and  Lumley''’  and  Huang.  Bradshaw  and  Coakley®).  Of  course,  for  non-equilibrium 
compressible  flows  involving  shocks  -  or  extremely  high  Mach  numbers  in  or  near  the  hyper¬ 
sonic  flow  regime  issues  related  to  near  wall  turbulence  modeling  are  even  more  unsettled 
(c.f.  Bnshnell' ). 

In  this  paper,  a  new  near-wall  two-equation  turbulence  model  of  the  K  —  £  type  is 
developed  for  high-speed  compressible  flows.  Two  features  distinguish  this  model  from  earlier 
work:  (a)  a  variable  density  extension  of  the  asymptotically  consistent  near-wall  model  of 
So  et  al.^  for  incompressible  flows  is  used,  and  (b)  high-Revnolds  nnmlx'r  models  for  tlu' 
dilatational  terms  are  implemented  as  developed  r<'C('ntly  by  .Sarkar  ('t  ah'  "*  and  Coleman 
and  Mansour"  based  on  an  analysis  of  homogeneous  turbnl('i\e('.  1  lu*  resulting  compp'ssibic 
model  is  test»'d  in  supersonic  flat  plate  boundary  hmus  at  zero  pressure'  gradient  both 


'odPS 


■s./or 


□  □ 


with  an  adiabatic  wall  and  with  wall  cooling  -  for  external  Mach  numbers  in  the  range 
0  <  Moo  <  10.  Comparisons  are  made  with  experimental  data,  and  with  the  predictions  of 
the  K  —  u  model  of  Wilcox^^,  in  an  attempt  to  address  the  following  questions: 

(a)  Can  incompressible  two-equation  models  be  extended  to  supersonic  boundary  layers 
with  a  minimum  of  compressible  corrections  that  are  systematically  determined? 

(b)  Will  these  models  perform  well  at  high  Mach  numbers? 

(c)  Where  will  the  models  break  down,  and  what  deficiencies  give  rise  to  it? 

These  issues  will  be  addressed  in  detail  in  the  sections  to  follow  and  recommendations 
will  be  made  for  future  research. 


where 

P  =  pRT  (4) 

2 

+  Uj,.)  (5) 

are,  respectively,  the  thermodynamic  pressure  and  viscous  stress  tensor  given  that  p  is  the 
density,  u,  is  the  velocity,  T  is  the  temperature,  R  is  the  ideal  gas  constant,  p  is  the  dynamic 
viscosity,  k  is  the  thermal  conductivity,  and  is  the  specific  heat  at  constant  pressure.  Here, 
the  Kinstein  summation  convention  applies  to  repeated  indices  and  (■),,  denotes  a  gradient 
with  respect  to  the  spatial  coordinate  x,. 
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Any  flow  variable  T  can  be  decomposed  into  mean  and  fluctuating  parts  in  two  ways: 

^  =  F  +  r,  T  =  (6) 

where  J-  represents  the  traditional  ensemble  average  whereas  J-  represents  a  Favre  average 
defined  by  _ 

The  mean  equations  of  motion  for  compressible  turbulence  take  the  form: 

+  ipui),i  =  0  (8) 

Q  _  _  _2_  _ 

—  (pUi)  +  =  -p,i  -  +  [p(n,,j  +  (9) 

d  d — 

—ipCpf)  +  {pu,Cpf),,  =  -^  +  “.p,.  +  u'iP^  +  “..j 


+cr,ju"j  +  pe  —  {pCpQ,),i  +  (/rT,),, 


where 


A,  =  u"u" 


Q,  =  u':T" 


=  --^puk.kStj  +  p{ui,j  +  iij.i)  (13) 

p€  =  a',y,j  (14) 

are,  respectively,  the  Reynolds  stress  tensor,  the  Reynolds  heat  flux,  the  mean  viscous  stress 
tensor,  and  the  turbulent  dissipation  rate.  Eqs.  (8)  -  (10)  are  derived  subject  to  one  major 
assumption:  turbulent  fluctuations  in  the  viscosity,  thermal  conductivity  and  specific  heat 
can  be  neglected. 

The  pressure  gradient-velocity  correlation  u'p(,  can  be  written  in  the  equivalent  form  (see 
Speziale  and  Sarkar*^) 

=  -(pRT<),,  -k  ipRuff"),,  -7<,-  (15) 

Consequently,  in  order  to  achieve  closure,  models  are  needed  for:  (a)  the  Reynolds  stress  Tij, 
(b)  the  Reynolds  heat  flux  Q,,  (c)  the  turbulent  dissipation  rate  t,  (d)  the  pressure  dilatation 
correlation  p'u',,  and  (e)  the  mass  flux  u".  Consistent  with  the  recent  work  of  Sarkar  et  al.®, 
the  dissipation  is  decomposed  into  solenoidal  and  compressible  parts: 


£  =  e.  +  e, 
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where  for  homogeneous  turbulence  'pe,  —  puj'-uj'-  and  pSc  —  given  that  tu'  is  the 

fluctuating  vorticity.  Here,  Ss  represents  the  dissipation  associated  with  the  energy  cascade. 
The  length  and  time  scales  will  be  built  up  from  £,  which  is  to  be  obtained  from  a  modeled 
transport  equation.  The  Sarkar  et  al.^  model  for  the  compressible  dissipation  will  be  used 
in  the  form: 

£c  =  otMfss  (17) 

where  Mt  =  {u" u'-  / ^ RT)^  is  the  turbulence  Mach  number  given  that  7  =  C'p/C'„  is  the 
ratio  of  specific  heats.  Here,  the  constant  a  is  approximately  0.5  based  on  direct  numerical 
simulations  of  homogeneous  shear  flow*®. 

The  Reynolds  stress  tensor  is  modeled  in  the  standard  eddy  viscosity  form 

nj  =  \l<i„  -  -iCJ.—iS,,  -  (18) 


where  K  =  It,,  is  the  turbulent  kinetic  energy,  5,j  =  ^{duifdxj  +  duj(dxi)  is  the  Favre- 
averaged  rate  of  strain  tensor,  Cp  is  a  dimensionless  constant  taken  to  be  0.096,  and  is 
a  wall  damping  function.  The  wall  damping  function  -  which  goes  to  one  sufficiently 
far  from  the  wall  -  will  be  discussed  in  more  detail  later.  A  simple  eddy  viscosity  model 
was  chosen  since  we  are  dealing  with  two-dimensional  attached  boundary  layers;  for  other 
applications  involving  more  complex  flows,  (18)  can  be  easily  generalized  to  an  algebraic 
stress  model  or  an  anisotropic  eddy  viscosity  model  (see  Rodi*'*  and  Speziale*^). 

The  Reynolds  heat  flux  and  mass  flux  terms  are  modeled  by  the  standard  gradient  trans¬ 
port  hypotheses: 


u'T" 


Ptt  e. 


(19) 


u. 


CJ,  K^_ 

- - P, 

P^p  ^  s 


(20) 


where  Prj  is  the  turbulent  Prandtl  number  which  is  taken  to  be  0.9  and  CTp  is  the  mass  flux 
constant  which  assumes  a  value  of  0.5.  A  recently  derived  model  by  Sarkar*®  is  used  for  the 
pressure  dilatation  correlation  which  is  given  by 

p'u[,  -  -axfPM^  -f  Q2pesMf  (21) 

where  V  =  is  the  turbulence  production.  Based  on  direct  numerical  simulations  of 

homogeneous  shear  flow,  Sarkar*®  determined  that  f>i  ~  0.4  and  Q2  «  0.2. 

In  order  to  achieve  closure,  modeled  transport  e(iuations  for  the  turbulent  kinetic  energy 
A’  and  the  solenoidal  part  of  the  turbulent  dissipation  rate  e*  are  needed.  The  exact  transport 
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equation  for  the  turbulent  kinetic  energy  is  given 


where 


+  {P^iK),i  =  -pTijUi,]  -  pe  +  pfu'i  i  -  u'lpi  +  u'Vij, j  +  Vli  + 


are,  respectively,  the  turbulent  and  molecular  diffusion  terms.  If  turbulent  fluctuations  in 
the  viscosity  are  neglected,  as  well  as  other  higher-order  terms,  the  molecular  diffusion  can 
be  approximated  as 

=  FA-,,  (24) 

in  the  boundary  layer.  Formally,  (24)  is  the  leading  order  part  of  (23)  close  to  the  wall 
provided  that  turbulent  fluctuations  in  the  viscosity  and  density  can  be  neglected.  The 
standard  gradient  transport  hypothesis  is  applied  to  the  modeling  of  the  turbulent  diffusion 


term; 


Vj  = 

<^K 


where  again  pr  =  pC^^f^K'^ le^  is  the  eddy  viscosity  and  is  a  constant  taken  to  be  0.75. 
This  leads  to  the  final  modeled  transport  equation  for  K : 

Q 

=  -pTijfx.j -p(l +aMf)£:5 

+Q2p£,Mf  -  +  (p  +  — )  K,i 

L  V  aA'  /  . 


The  exact  transport  equation  for  the  solenoidal  part  of  the  turbulent  dissipation  rate  is 
of  the  general  form: 

^(P£*)  +  (P“.£ J,.  =  -^pe^w.,,  +  P£3^^  T  Pc,  -  +  Pc,  +  (pe».,),i  (27) 

where  is  the  production  of  solenoidal  dissipation,  $i,  is  the  corresponding  destruction 
term,  and  is  the  turbulent  diffusion  of  solenoidal  dissipation.  For  the  sake  of  brevity,  the 
details  of  the  higher-order  correlations  that  comprise  Pf,,<^c,  ^-nd  Pj,  are  not  given  (these 
correlations  are  quite  complicated).  In  deriving  the  viscous  term  in  (27),  turbulent  fluctua¬ 
tions  in  the  molecular  viscosity  have  been  neglected  consistent  with  the  earlier  derivations. 
The  term  Pj,  in  (27)  represents  the  production  of  solenoidal  dissipation  by  deviatoric  mean 
strains  as  well  as  by  density  fluctuations.  This  term  will  be  modeled  as  a  variable  density 
extension  of  its  incompressible  form  with  density  fluctuations  neglected'^.  Hence,  we  take 

p.,  =  -CelP^T,j{ll,  j  -  (28) 
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where  C^i  is  a  dimensionless  constant  taken  to  be  1.5.  The  same  approach  is  used  in  the 
modeling  of  the  destruction  term  which  yields 

(29) 

where  /2  is  a  wall  damping  function  and  0^2  is  a  dimensionless  constant  which  is  taken  to 
be  1.83.  The  standard  gradient  transport  hypothesis  is  used  for  the  turbulent  diffusion  term 

!>..  =  _  (30) 

where  is  a  dimensionless  constant  which  is  taken  to  be  1.45. 

Modeling  the  variable  viscosity  term  in  (27)  is  a  bit  trickier.  For  the  high  Mach  number 
flows  to  be  considered  in  this  paper,  variations  of  the  viscosity  with  temperature  must  be 
accounted  for.  In  this  vain,  we  use  Sutherland’s  law*®  for  which 


( ^  ^  x/^o\ 
\i  +  x/f) 


(31) 


where  To  and  JIq  are  the  reference  temperature  and  viscosity  whereas  x  's  a  constant  that 
depends  on  the  gas  (for  air,  x  =  110°  A').  Eq.  (31)  could  also  be  used  for  the  viscosity  trans¬ 
port  term  (ljV)DV/Dt  in  (27).  However,  this  would  couple  the  dissipation  rate  transport 
equation  to  the  continuity  and  energy  equations  at  the  highest  time  derivative  -  an  undesir¬ 
able  feature  that  can  cause  numerical  stiffness.  Hence,  following  the  work  of  Coleman  and 
Mansour**,  we  will  use  the  power  law  approximation 


(where  n  ^  0.7)  for  the  formulation  of  the  viscosity  transport  term  {1(17)017/01  in  (27); 
Sutherland’s  law  will,  however,  be  used  for  the  calculation  of  the  viscous  diffusion  teiins  since 
it  is  more  accurate.  Coleman  and  Mansour**  showed  that,  for  an  isentropic  compression, 
(32)  yields  the  relation 

=  ^  =  [1  -  ”(7  -  !)]«.,.  (33) 

j/  Ot 

where  7  =  Cp/C\.  is  the  ratio  of  specific  heats.  On  physical  grounds  one  would  expect  (33) 
to  constitute  the  leading  order  part  of  {\ (17)017 ( Ot  and  it  will  be  used  to  avoid  numerical 
stiffness  problems.  The  final  form  of  the  modeled  dissipation  rate  equation  is  obtained  by 
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substituting  (28)  -  (33)  into  (27)  which  yields 


dt 


(/?£,)  +  {puies),i 


-  +  n(7  -  1) 


pE^U^^i  CelP  (ht,j 


(34) 


where  n  ^  0.7.  The  mean  viscosity  p  is  approximated  by  Sutherland’s  law  in  (34)  as  well  as 
in  the  other  transport  equations. 

3.  NEAR  WALL  MODIFICATIONS  FOR  COMPRESSIBLE  FLOWS 

Now  we  will  address  the  issue  of  near-wall  modeling  in  more  detail.  Two  fundamental 
assumptions  will  be  made:  (i)  the  functions  are  analytic  so  that  they  can  be  expanded  in  a 
Taylor  series  near  the  wall,  and  (ii)  turbulent  fluctuations  of  the  density  and  temperature 
vanish  at  the  wall.  While  the  former  assumption  is  well  accepted,  the  latter  one  -  despite 
its  widespread  use  -  is  somewhat  debatable.  However,  since  this  assumption  is  rigorously 
valid  for  isothermal  wall  conditions  -  and  a  good  approximation  for  the  adiabatic  wall  -  we 
are  justified  in  using  this  simplification  for  the  present  paper.  Without  this  assumption,  the 
near  wall  asymptotics  become  much  more  uncertain. 

A  Taylor  expansion  then  yields: 


p  =  diy  +  +  •  •  • 

(35) 

T"  =  'ny  +  72«/^  +  •  •  • 

(36) 

where  j3,  and  7,  are  functions  of  x.  ^  and  t  given  that  y  is  the  coordinate  normal  to  the 
wall.  Due  to  (35)  and  the  continuity  equation  (1),  it  follows  that  the  fluctuating  velocity 
components  (u,  i\w)  have  the  expansions 

u  —  a\y  -|-  Qiy^  -f  •  •  • 

(37) 

V  =  b^y'^  -1-  b-^y^  +  •  •  • 

(38) 

w  =  ciy  +  ciy'^  -!-••• 

(39) 

near  the  wall  based  on  either  standard  or  Favre  averages.  Eciuations 
to  their  incompressible  counterparts  a  simplification  that  arises 
vanishing  density  fluctuations  at  the  wall.  By  using  (37)  -  (39).  it  is 

(37)  -  (39)  are  identical 
from  the  assumption  of 
straightforward  to  show 

that  we  have  the  following  near- wall  asymptotics  for  the  crucial  turbulence  correlations: 

£,  =  0(1). 


K  =  0{y\ 


(40) 


yT  =  0{y%  p'<.  =  0(j/2) 

(41) 

{ufT"),  =  e)(y2),  (<),i  =  0{y^) 

(42) 

V,,=Oiy),  ^e.  =  0{l) 

(43) 

(the  precise  asymptotic  behavior  of  the  turbulent  diffusion  terms  (25)  and  (30)  are  not 
of  consequence  since  they  constitute  higher-order  terms  near  the  wall).  The  asymptotic 
constraints  (40)  -  (43)  can  be  satisfied  identically  with  only  two  damping  functions,  and 
/2,  which  behave  as 

L  =  h  =  0{y^)  (44) 


near  the  wall. 

Consistent  with  the  underlying  assumption  that  the  density  fluctuations  vanish  at  the 
wall  -  which  renders  the  fluctuating  velocity  to  be  solenoidal  to  the  leading  order  near  the 
wall  -  a  variable  density  extension  of  the  incompressible  wall  damping  functions  of  So  et 
al.®  wilt  be  used.  Validations  of  these  wall  damping  functions  have  been  carried  out  for  a 
number  of  two-dimensional  flows,  including  internal  and  external  flows  as  well  as  flows  with 
heat  transfer^^"^®.  In  all  of  these  cases,  good  correlations  with  direct  simulation  data^^"^'* 
and  measurements^®  have  been  obtained.  These  wall  damping  functions  are  implemented  in 
the  form: 


/m  =  tanh(y+/115) 

^3  L  ^€2  \^3^s/ 


(45) 

(46) 


where 


y'^  =  yUrli'  (47) 

U  =  R^  =  K^Ve,  (48) 

e:  =  -  2vh'/y^  (49) 

=  (50) 


In  (47)  -  (50),  1/  =  yip  is  the  local  mean  kinematic  viscosity  and  Ur  is  the  friction  velocity 
defined  by 

=  sJtuJK-  (51) 

where  is  the  wall  shear  stress^®. 
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4.  THE  FLAT-PLATE  BOUNDARY  LAYER  EQUATIONS 


The  steady  boundary  layer  form  of  the  mean  turbulence  equations  corresponding  to 
this  model  are  now  provided.  Consistent  with  most  practical  computations  of  compressible 
flows^^,  we  will  solve  the  energy  equation  in  its  total  enthalpy  form  where 


H  =  Cpf  +  -UiUi  +  K 


is  the  Favre-averaged  total  enthalpy.  The  resulting  boundary  layer  equations  at  zero  pressure 
gradient  take  the  form: 

^{pu)  +  ^{pv)  =  0  (53) 


_ du  _ du  d  \.  ^dii 


_,dH  _dH  d 
dx  dy  ~  dy 


J _ ]_\1  ^ 

<7 A-  Prj/J  dy 


_,dl<  _,dK 

"“to  + 


pu—  pv—  =  HT  -  /3(1  +  oiM;)es  -  otiPT  Mt 


+a,pe,M^+^i-\p^]+^  [Jl+^ 


d  (  du\  d 


dy  \  dy)  dy 


PT  \  dK 


(TkJ  d  y 


-,de,  _,de. 


1  ,  ,  dv\ 

-[3+n(l-l)],(^^  +  ^j.. 


^  e,  ( du 


^  ^  de. 


where  pj  =  pCpfpK^fCa  is  the  turbulent  viscosity.  A  summary  of  the  values  of  the  constants 
is  given  as  follows:  =  0.09^,  Cii  =  1.5,  Cj2  =  1-83,  =  0.75,  <Tj  =  1.45,  cTp  =  0.5,  a  = 

0.5,  oi  =  0.4,  Q2  =  0.2,  and  n  =  0.7.  The  damping  functions  are  as  specified  in  (45)  -  (46); 
consequently,  the  two  equation  model  of  So  et  al.®  is  recovered  in  the  incompressible  limit. 
The  turbulent  Prandtl  number  Prj  is  specified  as  0.9;  the  molecular  Prandtl  number  Pr  is 
taken  to  be  0.74  for  air  and  0.70  for  helium;  and  the  specific  heat  ratio  7  is  1.4  for  air  and 
1.67  for  helium.  Sutherland's  law  (31)  is  used  to  evaluate  the  molecular  viscosity  in  (54)  - 
(57)  for  air;  for  helium  the  power  law  quoted  in  Fernholz  and  F'inley^^  is  used. 


9 


A  few  comments  are  in  order  concerning  the  derivation  of  the  transport  equation  for  the 
total  enthalpy  H.  VVe  obtain  this  equation  by  adding  the  modeled  transport  equations  for 
6’pT,  and  K.  As  an  alternative  approach,  the  exact  transport  equation  for  //  can  be 

modeled  directly.  However,  we  do  not  feel  that  this  is  a  good  approach  since  //  is  not  a 
Galilean  invariant  quantity. 

Equations  (53)  -  (57)  are  solved  subject  to  the  standard  boundary  conditions: 

u  =  0,  V  =  0,  K  =  0,  ts  =  2'Uw{d\fK ! dy)^  (58) 

at  the  wall.  Two  different  types  of  thermal  wall  boundary  conditions  will  be  considered;  the 
adiabatic  wall  where  dll/Oy  is  zero  and  the  cooled  wall  case  where  H  is  constant.  .At  the 
edge  of  the  boundary  layer,  the  mean  velocity  and  the  total  enthalpy  are  required  to  match 
the  specified  free-stream  conditions.  On  the  other  hand,  the  turbulence  (piantities,  K  and 

are  assumed  to  be  zero  in  the  free  stream.  Thus  formulated,  the  above  equations  with  the 
appropriate  boundary  conditions  can  be  solved  numerically  using  the  boundary-layer  code 
developed  by  .Anderson  and  Lewis^'  and  modified  by  So  et  al.^ 

5.  DISCUSSION  OF  RESULTS 

Two  versions  of  tlu'  A'  — £  model  derived  herein  are  used  to  carry  out  the  compressible  fiat 
plate  boundary-layer  calculations  to  be  presented  in  this  section.  One  version  consists  of  the 
full  compr('ssible  model  as  given  in  Eqs.  (53)  -  (57)  and  is  designated  as  K  — £  model/1.  The 
si’cond  version  which  is  designated  as  A'-£  model/2  -  is  a  variable  density  extension  of  the 
incompressi!)le  two-equation  model  of  Reference  8  wherein  the  explicit  compressible  terms  in 
(56)  and  (57)  ar('  neglected  along  with  the  OKjOy  term  in  (55).  Furthermore,  comparisons 
will  also  be  made  with  results  of  the  well-known  l\  — -x!  model  of  Wilcox'^.  The  calculations 
cjf  these  two  different  versions  of  the  new  A  —  £  model,  as  well  as  those  for  the  A  —  ^ 
model,  will  allow  us  to  evaluate  the  range  of  validity  of  Morkovin's  hypothesis  and  assess 
the  im{)ortance  of  having  an  asymptotically  consistent  near-wall  correction  for  turbulence 
models.  Iln'iefore.  the  central  questions  po.sed  in  the  Introduction  can  be  addres.sed  by 
comparing  these  results  with  well-documented  <‘xp<'rimental  data  spanning  a  wide  range  of 
free-st if'am  Mach  numix'rs  and  cooled  wall  conditions. 

The  thif'c  diffi'rent  models  mentioned  above  are  u.sed  to  calculate  compressible  flat  plate 
bcaindary  layers  un  adia!)atic  as  well  as  cooled  walls  and  the  results  are  compared  with 
measurements  drawn  from  References  28  and  2f).  I  he  calculations  are  carried  out  over  the 
Mach  nuinl)er  range  0  <  <  10  for  the  adiabatic  wall  boundary  coiulition  and  over  the 

tem[)erature  range.  0  <  /  ,  //,,„  <  1  for  the  cooled  wall  case.  Here,  is  the  free-stream 
Mach  number.  I  ^  is  the  tr('<'-st r<'am  temperat ur(' and  is  the  adiabatic  wall  temperature. 
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Since  only  mean  flow  properties  are  available  from  Reference  28,  comparisons  are  made  with 
these  measured  quantities  only  and  the  reported  skin  friction  coefficient,  C/  =  2Tyj/p^U'^ 
(where  Uoo  is  the  free-stream  velocity  and  'p^  is  the  free-stream  mean  density).  In  addition, 
comparisons  of  the  calculated  turbulence  quantities  are  made  with  the  results  of  the  K  —  ui 
model.  All  calculations  are  carried  out  to  the  same  momentum  thickness  Reynolds  number 
{Re)  as  the  measurements  and  the  comparisons  are  made  at  these  respective  locations.  Four 
sets  of  data  are  chosen  from  Reference  28.  These  are  cases  55010504,  53011302,  and  73050504 
with  the  adiabatic  wall  boundary  condition  and  case  59020105  with  a  constant  cooled  wall 
temperature.  M^o  for  these  cases  are  2.244,  4.544,  10.31,  and  5.29,  respectively,  and  the 
corresponding  values  of  Re  are:  20,797;  5,532;  15,074;  and  3,939.  The  variations  of  Cj  with 
Moo  and  Too/Taw  are  compared  with  che  van  Driest  II  results  given  in  Reference  29.  Finally, 
the  mean  velocity  profiles  are  compared  with  the  van  Driest  law-of-the-wall  for  compressible 
boundary  layers^°’^‘  and  an  assessment  of  the  effects  of  compressibility  on  boundary-layer 
flows  is  attempted. 

The  mean  velocities  are  presented  in  two  dififerent  forms:  one  in  terms  of  wall  variables, 
u'^  versus  In  y'^,  and  another  in  terms  of  outer  variables,  u/Ux,  versus  y/6.  Here,  = 
^/^T)  and  6  is  the  boundary-layer  thickness.  In  this  way,  the  effects  of  the 

calculated  Ur  and  p  on  the  mean  velocities  can  be  assessed.  The  mean  temperature  is 
plotted  in  the  form  of  T /Tx  versus  y/S  while  all  turbulence  quantities,  such  as  K  and  — uu, 
are  normalized  by  uj  to  give  k'^  and  — uF*',  respectively,  for  presentation.  A  comparison 
with  the  van  Driest  law-of-the-walP°’^^  for  compressible  flows  is  also  attempted;  therefore, 
the  mean  velocities  of  two  cases  are  plotted  in  terms  of  the  compressible  defined  according 
to  Reference  31  as  =  /q^  du'^ .  Since  the  K  budgets  for  all  cases  considered  are 

essentially  similar,  only  the  budget  of  K  in  the  near- wall  region  for  case  73050504  (A/oo  = 
10.31  and  T^jTr  =  1.0)  is  presented.  Here,  the  reference  temperature  Tr  is  taken  to  be  the 
recovery  temperature  for  the  cases  with  an  adiabatic  wall  boundary  condition  and  taken  to 
be  the  adiabatic  wall  temperature  for  the  cooled-wall  case. 

The  u'^  versus  In  y^  plots  for  the  four  cases  are  shown  in  Figures  1  and  2.  In  each  figure, 
the  calculated  and  measured  skin  friction  coefficient  C/,  as  well  as  that  determined  from  the 
van  Driest  II  formula  of  Reference  29,  are  listed  for  comparison.  A  standard  law-of-the-wall 
given  by  u'^  =  (l/«:)/n  y^  +  B  \s  also  shown  for  comparison,  where  the  von  Karman  constant 
K  =  0.41  ha.s  been  assumed.  It  is  recognized  that  the  intercept  B  should  be  a  function  of 
A/oo,  however,  in  these  figures,  B  is  taken  to  be  4.7.  The  actual  value  used  here  is  not  too 
important  because  the  objective  is  to  determine  and  compare  the  variation,  if  any,  of  the 
log-law  slope  with  Mx  and  TuijTr.  The  log-law  is  seen  to  hold  true  for  all  cases  calculated 
using  both  the  K  ~e  and  K  —  u  models  (Figures  1  and  2)  and  the  value  of  k  thus  determined 
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for  the  three  adiabatic  wall  cases  is  essentially  0.41.  For  the  cooled- wall  case,  the  calculated 
K  is  not  equal  to  0.41  and  varies  from  model  to  model.  The  best  agreement  with  data  is 
given  by  the  calculation  using  K  —  e  model/2  while  the  smallest  value  of  «  is  predicted  by 
the  K  —u  model.  All  models  considered  yield  calculations  of  C/  that  are  in  agreement  with 
measured  data  and  with  van  Driest  II  values.  The  maximum  error  is  less  than  4%  (Figures  1 
and  2).  An  exception  is  the  cooled- wall  case  where  the  measurement  is  higher  than  the  van 
Driest  II  value.  .According  to  Reference  28,  the  measured  C/  in  this  case  is  not  as  accurate 
as  in  the  other  cases  and  this  explains  the  discrepancy  between  the  calculated  and  measured 
Cj  shown  in  Figure  2.  This  comparison,  therefore,  shows  that  the  additional  compressible 
terms  in  the  governing  equations  do  not  significantly  affect  the  calculated  results  so  long 
as  the  near-wall  model  used  is  the  same  and  is  asymptotically  consistent.  In  other  words, 
Morkovin’s  hypothesis  is  valid  for  flows  with  free-stream  Mach  numbers  as  high  as  10  and 
wall  temperature  ratios  noticeably  smaller  than  one. 

In  the  past,  velocity  profiles  in  wall  coordinates  were  invariably  plotted  in  terms  of 
to  illustrate  the  existence  of  the  van  Driest  log-law  and  the  constancy  of  k  in  compressible 
boundary-layer  flows®’^°’^k  Since  then,  the  compressible  law  of  the  wall  is  typically  taken 
to  be  given  by  u/"  rather  than  by  and  k  is  considered  to  be  about  0.41  and  constant 
over  the  Mach  number  range  of  0  <  A/oo  <  5.  The  calculated  and  measured  velocity  plots 
given  in  Figures  1  and  2  show  support  for  the  compressible  law-of-the-wall  when  it  is  written 
in  terms  of  u'^  rather  than  Furthermore,  k  is  determined  to  be  approximately  0.41 
and  is  relatively  constant  over  the  Mach  number  range  of  0  -  10  for  the  adiabatic  wall 
boundary  condition.  These  results  seem  to  conflict  with  the  proposal  of  van  Driest^°.  In 
order  to  resolve  this  seeming  contradiction,  the  velocity  plots  of  versus  In  for  cases 
5.5010504  (A/oo  =  2.244, Tu,/Tr  =  1)  and  .53011302  {Moo  =  4.544, Tu,/Tr  =  1)  are  shown  in 
Figure  3.  In  addition,  the  compressible  law-of-the-wall  as  given  in  Reference  31  is  shown  for 
comparison.  It  can  be  seen  that  a  line  that  is  parallel  to  the  compressible  law-of-the-wall 
can  be  drawn  through  a  few  of  the  date  points  spanning  over  a  narrow  range  of  y^.  On  the 
other  hand,  the  calculated  profiles  are  in  agreement  with  data  over  a  wider  range  of  yjJ;.  The 
slopes  of  the  calculated  profiles  are  roughly  parallel  to  that  determined  from  measurements 
and  are  slightly  larger  than  the  slope  of  the  compressible  law-of-the-wall  shown.  Therefore, 
irrespective  of  how  the  velocity  profiles  are  plotted,  the  calculations  are  in  good  agreement 
with  data.  However,  the  slope  of  the  log-law  appears  to  be  given  by  (0.41)-'  only  when  the 
profiles  are  plotted  in  terms  of  u'^ . 

Based  on  the  above  comparisons,  it  seems  that  there  is  very  little  difference  between 
the  predictions  of  the  K  —  e  and  K  —  u;  models.  This  is  particularly  true  in  regard  to  the 
calculations  of  Cj.  The  semi-log  plots  shown  in  Figure  1  tend  to  mask  the  differences  found 
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between  the  models  in  the  calculations  of  the  mean  temperature  and  density  profiles.  These 
differences  begin  to  show  up  in  the  plots  shown  in  Figure  3  and  in  the  predictions  of  the 
cooled-wall  case  (Figure  2).  In  order  to  determine  the  actual  difference  between  the  mean 
flow  predictions  of  the  K  —  e  and  K  —  u  models,  plots  of  ujUoo  versus  yjb  are  given  in  Figures 
4  and  5  for  the  adiabatic  and  cooled  wall  conditions,  respectively,  while  the  corresponding 
mean  temperature  profiles  are  plotted  in  Figures  6  and  7.  These  results  show  that  there 
are  substantial  discrepancies  between  measurements  and  the  profiles  of  ii/Uco  and  T /Too 
calculated  using  the  K  —  u  model.  The  discrepancies  increase  as  Moo  increases  and  as  T^/Tt 
decreases.  The  best  results  are  given  by  A'  —  £  model/1.  On  the  other  hand,  the  predictions 
of  A'  —  £  model/2  are  very  close  to  those  of  K  —  e  model/ 1  and  are  substantially  different 
from  those  of  the  K  —  ui  model.  This  is  further  evidence  that  an  asymptotically  consistent 
near-wall  model  is  more  important  than  the  additional  compressible  terms  as  far  as  the 
predictions  of  the  mean  flow  properties  are  concerned. 

The  ability  of  the  K  —  ui  model  to  predict  the  variations  of  C/  with  Moo  and  T^/Taw  is 
well  established^^.  If  the  proposed  K  —  e  model  is  to  be  accepted,  its  ability  to  predict  Cj 
for  different  Moo  and  T^/Taw  has  to  be  demonstrated.  In  Figure  8,  the  variation  of  Cf/{Cj)i 
with  Moo  for  the  case  of  the  adiabatic  wall  boundary  condition  is  shown.  Here,  (C/),  is  the 
skin  friction  coefficient  for  an  incompressible  flow  evaluated  at  Re  =  10“*  which  is  determined 
to  be  2.70  X  10“^.  This  figure  shows  a  comparison  of  the  calculations  of  K  —  e  model/1  and 
K  —  £  model/2  with  the  van  Driest  II  curve''^^.  Essentially,  there  is  no  difference  between  the 
predictions  of  K  —  e  model/1  and  A  —  £  model/2  and  both  results  are  in  excellent  agreement 
with  data.  The  predictions  for  the  cooled  wall  case  are  shown  in  Figure  9.  Calculations  for 
this  case  are  carried  out  for  Moo  =  5  and  Re  =  10“*  where  the  incompressible  C j  is  again 
determined  to  be  2.70  x  10“^.  Three  sets  of  calculations  are  presented.  These  are  for  K  —  £ 
model/1,  A— £  model/2  and  a  third  version  of  K —e  model/1  where  the  dK/dy  term  in  (56)  is 
neglected.  It  can  be  seen  that  an  error  of  5%  or  larger  starts  to  accumulate  at  approximately 
Txu/Taw  =  0.4  for  K  —  e  model/1.  This  trend  is  contrary  to  the  results  reported  in  Reference 
12.  An  examination  of  the  governing  equations  solved  by  other  researchers  revealed  that, 
besides  the  differences  noted  in  the  turbulence  model  equations,  the  mean  enthalpy  equation 
solved  by  these  researchers  does  not  include  the  term  dK/dy  on  the  right  hand  side  of  (55). 
Indeed,  when  the  dK/dy  term  is  neglected,  an  overall  important  improvement  is  obtained. 
The  predicted  C/  at  Tu,/Taw  =  0.2  is  increased  by  about  6%,  thus  giving  a  much  better 
agreement  with  the  van  Driest  II  curve^^.  If  the  additional  compressible  terms  in  the  K  —  e 
equations  are  further  neglected  (K  —  £  model/2),  the  calculated  Cj  is  only  improved  by 
about  2%.  The  remaining  disagreement  could  he  attributed  to  the  assumption  of  a  constant 
turbulent  Prandtl  number  which  cannot  properly  account  for  a  reduction  of  turbulent  mixing 
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resulting  from  a  cooled  wall  (in  fact,  the  dK/dy  term  can  be  formally  dropped  if  Pr^  =  0.74; 
such  a  lowering  of  the  turbulent  Prandtl  number  would  not  be  inconsistent  with  experiments 
for  the  cooled  wall  case).  However,  this  effect  can  be  appropriately  incorporated  in  a  heat- 
flux  model  and  modifications  can  be  formulated  to  describe  its  near-wall  behavior^^’^^.  In 
other  words,  if  highly  cooled-wall  flows  are  to  be  predicted  correctly,  turbulent  heat  fluxes 
and  their  near-wall  behavior  need  to  be  modeled. 

The  distributions  of  k'^  and  —uv '  across  the  boundary  layers  are  compared  in  Figures 
10  -  13.  In  these  figures,  only  the  results  oi  K  —  e  model/1  and  the  K  —  u>  model  are 
compared.  The  results  for  the  adiabatic  wall  boundary  condition  are  shown  in  Figures  10 
and  12  while  those  for  the  cooled-wall  case  are  plotted  in  Figures  11  and  13.  It  can  be  seen 
that  the  K  —  u>  model  underpredicts  k'^  in  the  near-wall  region  and  overestimates  its  value 
in  the  outer  part  of  the  boundary  layer.  This  incorrect  prediction  is  common  for  all  cases 
studied.  The  predicted  near- wall  behavior  is  substantially  different  from  that  given  by  K  —  e 
model/1.  Instead  of  yielding  a  finite  slope  for  k'*'  at  the  wall,  a  near  zero  slope  is  calculated. 
Furthermore,  the  maximum  k'^  calculated  is  about  half  that  given  by  K  —  e  model/ 1.  On  the 
other  hand,  the  K  —  to  model  yields  the  correct  near-wall  behavior  for  —uiP'  but  overpredicts 
its  value  in  the  outer  part  of  the  boundary  layer.  The  overprediction  extends  across  the  range, 
0.2  <  y/6  <  1.0.  It  is  now  clear  that  the  K  —  u)  model  is  formulated  to  give  correct  results  for 
the  mean  velocity  and  the  wall  shear  stress;  however,  its  predictions  of  other  properties  are 
in  doubt,  particularly,  in  regard  to  the  near-wall  behavior.  Reduction  of  turbulence  activity 
in  the  outer  part  of  the  boundary  layer  is  clearly  evident  when  either  compressibility  or  wall 
cooling  effects  are  present.  The  reduction  increases  as  Moo  increases  and  decreases. 

This  is  further  substantiated  by  the  very  significant  drop  in  the  maximum  value  of  k'^  as  Moo 
increases.  Therefore,  it  is  expected  that  turbulence  activity  will  be  substantially  reduced  in 
a  flow  where  the  free-stream  Mach  number  is  large  and  the  wall  is  highly  cooled. 

The  predictions  of  the  near-wall  flow  can  be  further  examined.  With  the  assumption 
of  vanishing  p'  and  T"  at  the  wall,  Taylor  expansions  of  the  fluctuating  velocities,  density 
and  temperature  about  y  =  0  are  given  by  (35)  -  (39).  These  expansions  together  with 
the  definition  for  £,  can  be  rearranged  to  give  the  following  dimensionless  expansions  for 
k'^ ,  —uv  *  and  i.e., 


(59) 

(60) 

£+  =  2ak  +  46iy+  -t- . . . , 

(61) 

where  the  a’s  and  b's  are  time-averaged  coefficients  that  are  functions  of  x  and  From 
these  expansions,  it  can  be  easily  deduced  that  k^ / {y^Y  =  0.5  at  the  wall.  Therefore,  the 
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asymptotic  behavior  of  k'^ is  0.5  and  is  independent  of  Moa  as  well  as  the  thermal 
wall  boundary  conditions.  The  accuracy  to  which  a  model  can  predict  this  quantity  is  a 
reflection  of  its  asymptotic  consistency.  Finally,  the  variation  of  at  the  wall  with  Mach 
number  can  be  deduced  from  a  plot  of  Ufe  versus  Mool  therefore,  this  result  is  shown  in 
Figure  14  for  the  cases  with  the  adiabatic  wall  boundary  condition.  The  value  of  Qk  for  the 
corresponding  incompressible  flow  is  taken  from  Reference  8.  Model  calculations  of  ak,auv 
and  k'^ / are  tabulated  in  Table  1  for  comparison.  Since  the  K  —  jj  model  is  not  an 
asymptotically  consistent  near-wall  model,  its  predictions  of  these  limiting  values  are  poor; 
therefore,  they  are  not  listed  in  Table  1.  From  these  results,  it  can  be  seen  that  the  present 
K  —  £  model  is  asymptotically  consistent  up  to  a  Mach  number  Moo  of  10  for  both  adiabatic 
and  cooled  wall  conditions.  In  other  words,  the  predictions  of  k'^  and  — uir+  shown  in  Figures 
10  -  13  are  more  likely  to  be  correct  compared  to  those  given  by  the  A'  —  u!  model.  These 
results  further  show  that  ak  is  a  function  of  the  Mach  number  and  decreases  with  increasing 
Moo  (see  Figure  14).  This  means  that  viscous  dissipation  at  the  wall  also  decreases  with 
increasing  Moo- 

The  near-wall  budget  of  K  for  the  Moo  =  10.31  case  (73050504)  is  shown  in  Figure 
15.  It  is  evident  that  the  budget  of  K  bears  a  strong  resemblence  to  that  calculated  for 
incompressible  flows^^.  The  additional  compressible  terms  have  a  negligible  effect  on  the 
near-wall  K  budget.  Therefore,  the  assumptions  made  to  extend  the  near-wall  damping 
function  /2  in  the  dissipation-rate  equation  to  compressible  flows  appear  to  be  justified. 
Again,  it  can  be  seen  that  viscous  diffusion  balances  dissipation  at  the  wall.  This  balance 
extends  to  about  j/^  =  4  where  turbulent  diffusion  and  production  become  important.  In  the 
region,  4  <  <  15,  viscous  and  turbulent  diffusion  as  well  as  production  and  dissipation 

are  equally  important.  Beyond  y^  =  15,  production  and  dissipation  are  approximately  in 
balance,  just  as  in  the  case  of  incompressible  flows.  Consequently,  the  near- wall  behavior  of  K 
is  very  similar  for  both  incompressible  and  compressible  flows  (thus,  explicit  compressibility 
effects  are  not  that  important  in  the  near-wall  region). 

6.  CONCLUSIONS 

In  contrast  to  the  conclusions  drawn  by  Bradshaw  et  al.^  and  Huang  et  al.^,  the  present 
investigation  shows  that  conventional  K  —  e  models  can  be  extended  to  calculate  equilib¬ 
rium  compressible  boundary-layers,  at  high  Mach  numbers,  if  the  near-wall  modifications  to 
these  models  are  asymptotically  correct  and  internally  consistent.  When  properly  modified, 
K  —  £  type  models  can  be  used  to  calculate  equilibrium  compressible  boundary-layers  with 
free-stream  Mach  numbers  as  high  as  10  and  wall  temperature  ratios  as  low  as  0.2.  Two 
different  versions  of  the  K  —  £  model  proposed  in  this  study  are  tested:  one  with  all  of 
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the  systematically  derived  dilatational  terms  included  in  the  governing  equations  {K  —  e 
modcl/1)  and  another  which  is  a  variable  density  extension  of  the  incompressible  limit  of 
this  model  where  explicit  compressible  terms  are  neglected  (K  -  e  model/2).  Calculated 
mean  flow  properties  and  wall  shear  stresses  from  these  two  versions  of  the  K  —  e  model  are 
in  excellent  agreement  with  measurements.  The  near-wall  behavior  of  the  calculated  tur¬ 
bulence  properties  is  consistent  with  exact  results  and  the  correct  limiting  values  for  these 
properties  are  recovered.  Compressibility  effects  are  found  to  have  a  negative  influence  on 
turbulent  mixing  which  is  reflected  in  a  reduction  of  the  maximum  value  of  the  turbulent 
kinetic  energy  and  the  value  of  the  viscous  dissipation  at  the  wall  as  the  free-stream  Mach 
number  increases.  Wall  cooling  also  produces  the  same  effects  and  the  two  versions  of  the 
proposed  A  —  £  model  mimic  the  trend  very  well.  This  means  that  Morkovin’s  hypothesis 
is  valid  for  the  full  range  of  Mach  numbers  and  wall  temperatures  tested,  provided  that  the 
near-wall  two-equation  models  used  to  close  the  governing  equations  are  asymptotically  cor¬ 
rect  and  internally  consistent  in  the  near-wall  region.  On  the  other  hand,  the  well  accepted 
h  —  u)  model  provides  good  correlations  with  measurements  for  the  mean  velocity  and  wall 
shear  stress  only.  Its  predictions  of  the  temperature  profiles  are  in  considerable  error,  partic¬ 
ularly  for  high  free-stream  Mach  numbers  and  for  the  cooled-wall  boundary  condition.  The 
reason  for  these  discrepancies  is  traced  to  the  model’s  inability  to  reproduce  the  near-wall 
turbulence  properly;  the  modeled  asymptotic  behavior  of  the  turbulence  properties  is  not 
consistent  with  the  exact  equations  governing  the  transport  of  these  quantities.  Despite 
these  deficiencies,  it  is  still  rather  surprising  how  well  the  K  —  lo  model  predicts  the  mean 
velocity  and  skin  friction  at  high  Mach  numbers. 

Considering  the  excellent  results  obtained  in  this  paper  for  equilibrium  turbulent  bound¬ 
ary  layers,  future  applications  of  this  compressible  K  —  £  model  are  planned  for  the  study  of 
more  complex  non-equilibrium  boundary-layer  flows  involving  shocks.  In  these  applications, 
we  would  expect  the  explicit  compressible  terms  appearing  in  the  modeled  transport  equa¬ 
tions  for  K  and  e,  to  play  a  more  important  role.  After  these  further  tests  are  completed, 
we  will  have  a  much  better  idea  of  the  full  range  of  applicability  of  this  new  compressible 
K  —  £  model  for  wall-bounded  turbulent  flows. 
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Case  Moo  T-uj/Tr  Model 


55010504  2.244  1.0  A' -  e  model/1  0.0998  7.167 


1.0  A' -£  model/2  0.0992  7.198 


1.0  A' -£  model/1  0.0850  6.700 


53011302  4.544  1.0 


73050504  10.31  1.0  A' -£  model/1  0.0785  6.630 


73050504  10.31  1.0  A' -£  model/2  0.0771  6.740 

59020105  5.29  0.92  A' -e  model/1  0.0805  6.120 


59020105  5.29  0.92  A' -£  model/ 2  0.0788  6.140 


Cut,  X  10“* 

7.167 

0.50 

7.198 

0.50 

6.700 

0.50 

6.760 

0.50 

6.630 

0.50 

6.740 

0.50 

6.120 

0.50 

6.140 

1 

0.50 

Table  1.  Asymptotic  near-wall  behavior  of  the  calculated  turbulence  properties. 
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Figure  1.  Semi-log  plots  of  u"*"  for  adiabatic  wall  boundary  condition;  (a)  Moo  =  10.31,  (b) 
Moo  =  4.544,  and  (c)  A/oo  =  2.244. 
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Figure  3.  Semi-log  plots  of  ;/+  for  two  different  (a)  =  4.544  and  (b)  M 
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Figure  4.  Linear  plots  of  ujV.^  for  adiabatic  wall  boundary  condition:  (a)  =  10.31,  (b) 

=  4.541.  and  (c)  ~  2.244. 


Figure  6.  Linear  plots  of  T for  adiabatic  wall  boundary  condition:  (a)  Moa  —  10.31,  (b) 
Moo  =  4.544,  and  (c)  Moo  =  2.244. 
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Figure  7.  Linear  plots  of  T /Too  for  cooled  wall  boundary  condition. 
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Figure  9.  Variation  of  Cjl{Cj)i  with  for  A/^o  =  5.0. 
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Figure  10.  Behavior  of  k'^  across  the  boundary  layer  for  adiabatic  wall  boundary  condition 
(a)  A/oo  =  10.31,  (b)  A/oo  =  4.544,  and  (c)  A/oo  =  2.244. 
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Figure  11.  Behavior  of  k'^  across  the  boundary  layer  for  cooled  wall  boundary  condition. 


31 


Figure  12.  Behavior  of  -uF*"  across  the  boundary  layer  for  adiabatic  wall  boundary  condi 
tion:  (a)  =  10.31,  (b)  =  4.544,  and  (c)  Moo  =  2.244. 


34 


Loss  Gain 


35 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No  0704-0188 

PijDf'C  '•  ;  oofden  this  oi'e-TTion  .♦  oi.-,rfT>3i,  «.  ->si : oe'  ^espcrse  mciudihg  th^  tiTie  lOf  re-ziewng  ihstruaions  searching  e»'Stirig  data  vDufCPii 

gatherihd  jnd  -"airtainina  the  data  neeO'KJ.  ana  :  :.rpDierih  j  ’e'.>en,-ng  nlieaion  o*  tnior'r'alion  Send  comments  fe-aaroinq  this  burden  estimate  or  an»  othier  asoert  ot  this 

Olie't' jn  .1*  Tt  '.r  matron,  .n.,  luQing  suggi^stions  ^C'  reau^  ih ’h-s  Our-aen  t :  Washington  rieadQuarters  Ser#,ces  L"  rector  ate  ror  -nf  :>fmat'Cr'  Operations  and  Reports  U  1 S  .etterson 
DavisHijh/sa,  •suitei2C4  uflmgtor'  .i  «l  2?^02-4iQ2  and  to  mo  O't"  e  o»  N'ar^agemer't  and  dudget  Paperi^jOf*  ReOomon  Project  (0  704-0  ’  88)  Washir'gton,  DC  20S03 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3  REPORT  TYPE  AND  DATES  COVERED 

November  1991 _ Contractor  Report _ 

4.  TITLE  AND  SUBTITLE 

A  Near-Wall  Two-Equation  Model  for  Compressible 

Turbulent  Flows 

S.  FUNDING  NUMBERS 

NASl-18605 

505-90-52-01 

6.  AUTHOR(S) 

H.  S.  Zhang,  R.  M.  C.  So,  C.  G.  Speziale  and  Y.  G-  Lai 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Institute  for  Computer  Applications  in  Science 
and  Engineering 

Mail  Stop  132C,  NASA  Langley  Research  Center 

Hampton,  VA  23665-5225 

8.  PERFORMING  ORGANIZATION 

REPORT  number 

ICASE  Report  No.  91-82 

9.  SPONSORING  - MONITORING  AGENCY  NAME(S)  ANO  AOORESS(ES) 

National  Aeronautics  and  Space  Administration 

Langley  Research  Center 

Hampton,  VA  23665-5225 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

NASA  CR-189565 

ICASE  Report  No.  91-82 

11.  SUPPLEMENTARY  NOTES 

Langley  Technical  Monitor:  Michael  F.  Card  To  be  submitted  to  AIAA 

Final  Report  Journal 

12a  DISTRIBUTION  AVAILABILITY  STATEMENT 

Unclassified  -  Unlimited 

Subject  Category  34 

12b  DISTRIBUTION  CODE 

13.  abstract  (Maximum  200  words) 

A  near-wall  two-equation  turbulence  model  of  the  K-e:  type  is  developed  for  the  des¬ 
cription  of  high-speed  compressible  flows.  The  Favre-averaged  equations  of  motion 
are  solved  in  conjunction  with  modeled  transport  equations  for  the  turbulent  kinetic 
energy  and  solenoldal  dissipation  wherein  a  variable  density  extension  of  the  asymp¬ 
totically  consistent  near-wall  model  of  So  and  co-workers  is  supplemented  with  new 
dilatational  models.  The  resulting  compressible  two-equation  model  is  tested  in  the 
supersonic  flat  plate  boundary  layer  —  with  an  adiabatic  wall  and  with  wall  cooling 
—  for  Mach  numbers  as  large  as  10.  Direct  comparisons  of  the  predictions  of  the  nev 
model  with  raw  experimental  data  and  with  results  from  the  K-to  model  Indicate  that  It 
performs  well  for  a  wide  range  of  Mach  numbers.  The  surprising  finding  is  that  the 
Morkovin  hypothesis,  where  turbulent  dilatational  terms  are  neglected,  works  well  at 
high  Mach  numbers  provided  that  the  near  wall  model  is  asymptotically  consistent. 
Instances  where  the  model  predictions  deviate  from  'he  experiments  appear  to  be  at¬ 
tributable  to  the  assumption  of  constant  turbulent  1 randtl  number  —  a  deficiency 
that  will  be  addressed  in  a  future  paper. 

14  SUBJECT  terms 

compressible  turbulence  modeling;  near-wall  two-equation  models; 
supersonic  boundary  layers 

15  NUMBER  OF  PAGFS 

37 

16  PRICE  CODE 

AO  3 

17  SECURITY  classification 

OF  report 

Unclassified 

18  security  classification 

OF  THIS  PAGE 

Unclassified 

19  security  CLASSIFICATION 

OF  ABSTRACT 

20  LIMITATION  OF  ABSTRACT 

Std'^da^d  298  (9ev  2-89) 


Ov  'INS-  -ltd  MM- 'ft 


NASA-Langley,  1991 


7540-0'  280  S800 


